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0\ : ABSTRACT 

We study the effect of anisotropic diffusion on the one-dimensional annihilation 

O ; 

reaction kA— nnert with partial reaction probabilities when hard-core particles meet 
t^- ■ in groups of k nearest neighbors. Based on scaling arguments, mean field approaches 

and random walk considerations we argue that the spatial anisotropy introduces no 
^ . appreciable changes as compared to the isotropic case. Our conjectures are supported 

q ' by numerical simulations for slow reaction rates, for k = 2 and 4. 
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1. INTRODUCTION 

In this work we consider the effect of anisotropy on one-dimensional annihilation 
reaction processes kA— »inert, with diffusion and hard-core particle interactions. Recent 
interest in stochastic dynamics of low-dimensional many-body systems has been largely 
due to the fact that mean-field "rate equation" description applies in many instances 
not only in high dimensions, D, but also down to D — 1 or D = 2. Thus, in order 



to observe fluctuation-dominated behavior, one has to consider low-dimensional sys- 
terns. Specifically, in ID, a host of new exact solutions, asymptotically exact scaling 
arguments, and equivalence to other more traditional ID many-body systems, for in- 
stance, Heisenberg spin chains, have been reported [1-15]. Some of these results have 
experimental relevance [16-17]. 

Recent work by Janowsky [18] provided an indication that spatial anisotropy in 
diffusion rates (i.e., unequal hopping probabilities to the left and to the right on the 
ID line) modifies the fluctuation behavior of the annihilation reaction A+B— ^inert. In 
particular, the exponent of the large-time, t, power-law decay of the particle concentra- 
tion, 
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c(t) ~ t 



— a 



(1.1) 



which applies provided the initial A- and B-densities were equal, was changed from 
a = l/4tol/3. Janowsky's observation was largely numerical, and we are not aware of 
any analytical or phenomenological explanation available in the literature. 

In order to gain insight into the origins of the asymptotic behavior for anisotropic 
diffusion, we propose, in the present work, to consider reactions kA— nnert, i.e., k- 
particle annihilation on the line. The reason for favoring these reactions is that for 
k > 3 the asymptotic large-time behavior in the isotropic case is mean-field [3,19]. 
There are several phenomenological arguments for this mean-field behavior, as well as 
for the fluctuation-dominated behavior for k = 2 (while the case k = 3 is marginal). 
We survey these arguments in Section 2 and propose extensions to the anisotropic case. 

We then report, in Section 3, numerical results for k = 2 and 4. Both analytical 
considerations and numerical results suggest that spatial anisotropy has little effect on 
the overall fluctuation vs. mean-field regimes for these reactions. Large-scale numerical 
simulations indicate that details of the variation of the particle density with time are 
similar and allow verification of the exponent predictions. A brief summary is given is 
Section 4. 

2. PHENOMENOLOGICAL CONSIDERATIONS 

The early arguments [20] yielding the fluctuation-dominated, diffusive behavior, for 
instance, for the reaction 2A— nnert, were simplistic but powerful. Thus the line of the 
argument went something like this: the average particle density c(t) at time t implies 
the average separation of order 1/c (in ID). Given the time t, the diffusion constant 
(of a single particle) T>, and this separation, the only dimensionless combination is 
T>tc 2 , which implies c ~ (T>i)~ 1 / 2 . This approach was further refined in [3]. Indeed, 
assuming that the fluctuation behavior is universal, i.e., that the initial density p = c(0) 
is "forgotten" at large times, one can write down the large-time scaling relation (given 
here for ID only), 



where the scaling function F is 1 for small arguments but must be power-law for large 



c ~ pF{Vtp 2 ) , 



(2.1) 
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arguments such that the p-dependence is canceled out. The resulting prediction is 



These arguments are valid provided the initial configuration is random (has no 
correlations, so that the only length scale is 1/p). Their advantage is that they can be 
extended to D > 1 and, occasionally, to some more complicated reactions [3,21]. When 
compared to the predictions of the mean-field rate equations, they yield the upper 
critical dimension values. For instance for the k = 2 reaction considered earlier, one 
gets a = D/2 and, when compared to the mean-field value a = 1, the upper critical 
dimension is identified as D = 2. The disadvantage of such considerations is that they 
do not yield any associated approximation scheme for more detailed calculations. 

We will term the preceding arguments "the scaling approach." Making diffusion 
anisotropic modifies the diffusion constant and also introduces a new dimensional quan- 
tity, the drift velocity. Its effect on the scaling-approach predictions is not clear in 
general. We note, however, that the "scaling approach" was mainly used for the sim- 
plest two-body reactions for which the drift should not be important. Indeed, in the 
fluctuation-dominated regime the reaction rate eventually renormalizes to the "fast- 
reaction" (so-called "diffusion-limited") limit. The interparticle distributions are be- 
lieved to approach diffusion-dominated forms [13,22-24], with effects of other particle 
interactions (such as hard-core, etc.) being "renormalized away." In the drifting refer- 
ence frame, therefore, the reaction will be the same with slowed-down diffusion. Thus 
the results of [18] for the reaction A+B— >inert, which indicate that anisotropic diffu- 
sion can modify the fluctuation-dominated behavior, were surprising and they are still 
largely open to interpretation. 

Another approach was developed in [19]. Assuming that mean-field applies, one 
uses the uncorrelated form of the interparticle distance distribution to make specific 
predictions of the reaction event rates. The result is a mean-field rate equation. When 
combined with phenomenological arguments similar to those described earlier, i.e., dif- 
fusive time-scale estimates, this approach can predict, for instance, that the k > 3 
annihilation reactions are mean-field in ID (for isotropic diffusion), while the k = 3 
reactions are marginal, which was recently confirmed numerically and analytically [25], 
and the k = 2 case is fluctuation-dominated. Similar mean-field approximations for 



c(t) 



universal constant 



(2.2) 
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other ID correlation quantities were used to study more complicated mostly two-body 
reactions [23-24] in ID. The advantage of these approaches is that they are associated 
with specific calculational schemes which yield results in the mean-field approximation 
by typically providing a rate equation for the density and other quantities. 

The disadvantage of these mean-field based approximations is that they cannot be 
easily extended to D > 1 because explicit results for uncorrelated particle systems are 
used [19], or objects which are only appropriate in ID are considered [23-24]. (There 
are, however, numerous other mean-field formulations not limited to ID.) We note that 
the approach of [19] for instance, only requires that for fast diffusion (i.e., assuming 
negligibly slow reaction rate), the particles become uncorrelated. As in [19], we consider 
here particles diffusing on the line and reacting in groups of k on encounters. However, 
if less than k particles meet, they interact as hard-core objects. In Figure 1 we show a 
snapshot of the stochastic evolution of such an anisotropic process for k = 4. Further 
details will be given in Section 3. 

An important observation is that making the diffusion anisotropic does not change 
the property of the hard-core particle system to loose correlation at large times [4] , due to 
diffusion only (disregarding the reaction). This would suggest that the approximations 
of [19] should remain valid. The behavior for k > 3 should be mean-field for anisotropic 
hopping, with k = 3 marginal, while k = 2 is non-mean-field (although the arguments 
of [19] cannot predict the exponent a). 

Finally, we consider yet another line of argument based on more local, "relative 
coordinate" considerations. These ideas are similar to those advanced in the exact 
asymptotic-limit studies of two-particle reactions [11,13,22]. For illustration, let us 
consider the isotropic-case /c-particle annihilation in ID. We note that when k particles 
needed for reaction end up as a close group due to diffusion, the "memory" of this 
event will be "washed away" by diffusion provided the random walk in the k — 1 relative 
distances is non-recurrent. Only for the 1- or 2-dimensional relative-coordinate walk, 
will the encounter be likely repeated with the same particles (assuming that the system 
is already quite dilute). Thus, local fluctuations are likely to be less important for 
k — 1 > 2, which is consistent with the earlier identification of the borderline value 
k = 3 which separates the mean-field and fluctuation-dominated regimes. 

The advantage of the "relative coordinate" argument is that it is local and therefore 
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it can be easily applied in D > 1, etc. However, the "locality" is also the disadvantage. 
Indeed, the fluctuation-dominated behavior is generally considered to be a genuinely 
many-body effect. Few-particle arguments can at best provide hints on the various 
regimes of behavior. 

Since the hard-core constraint is short-range, it should have no effect when the 
particles are well separated, in the dilute limit. Therefore, the recurrence or nonrecur- 
rence property in the relative coordinates will be the same for isotropic or anisotropic 
diffusion (as long as no finite-size effects are present, i.e., as long as the particles cannot 
interact "around" the system for periodic boundary conditions, for instance). In fact, 
for k = 2 we anticipate small or no changes as compared to the isotropic case. Our nu- 
merical results (next section) have verified this expectation and also indicated a similar 
property for k = 4. 

3. NUMERICAL RESULTS 

The specific details of the hard-core interaction and annihilation reactions, used 
in our numerical simulations, are given below. We note, however, that the conclusions 
regarding the asymptotic large-time behavior are generally expected to be universal. 
The numerical results reported were "large-scale," computer-resource demanding, which 
was the main reason for considering only the cases k = 2 (non-mean-field) and k = 4. 
The latter provides the simplest example of the "clean" mean-field behavior already 
well-studied for the isotropic diffusion [19]. 

At time t = 0, each site of the one-dimensional lattice with periodic boundary 
conditions is occupied with probability c(0) = p by identical hard-core particles. The 
particles perform a biased random walk between nearest-neighbor lattice sites. At a 
given time t one of the N(t) particles present in the system is picked at random. Let 
j denote its lattice site location. This particle hops to the site j + 1 (or j — 1) with 
probability h (or 1 — h) provided the target site is vacant, otherwise it remains in 
place. After each successful hopping attempt, the active particle can annihilate with 
probability q < 1 with k — 1 consecutive particles located in the direction of the hopping, 
i.e., on sites j + 2, . . . , j + k (or j — 2, . . . ,j — k), provided of course that all the 
k — 1 corresponding neighboring sites were already occupied. If the target site j + 1 
(or j — 1) was already occupied, so that the hopping event did not take place, the 
active particle at j may still annihilate (with probability q) with k — 1 particles at sites 



- 5 - 



j + 1, . . . , j + k — 1 (or j — 1, ... ,j — k + 1), provided they were all occupied. This rule, 
involving annihilation with particles in the direction of the hopping attempt, successful 
or unsuccessful, introduces correlation between hopping and reaction which has some 
similarity with the actual chemical reactions in D > 1; see [19]. 

Each successful annihilation reduces the total number of particles present in the 
system by k. The numerical procedure allows for N such hopping with reaction attempts 
per each time unit. Thus, only after iV attempts the time is increased by 1, and the 
particle number N(t + 1) is recalculated. This methodology is particularly efficient for 
dilute regimes and sets up a well-defined time scale since on average each particle is 
selected once per unit time. Of course, the results are only meaningful for sufficiently 
large system sizes so that N(t) 3> k holds for all times t studied. (Strictly speaking, 
low reaction probability is also required when particle density is of order 1.) 

Our Monte Carlo simulations confirm the general theoretical expectations for k = 2 
and 4. Results for the density, c(t), are illustrated in Figure 2. In order to investigate any 
changes due to anisotropy of hopping in the large-time asymptotic regime, we studied 
the case h = 0.5, compared to the maximal anisotropy, h = 1. In fact, large-scale 
simulations up to 10 8 time steps, starting with random, homogeneous initial particle 
distribution with densities p = 0.25, 0.5 and 0.8, indicate that no such differences are 
present, even for rather slow reactions, g -C 1 (our numerical values were as low as 
0.001). Moreover, the evolution of the particle concentration for h = 1 closely follows 
the density of the isotropic case for all times. For instantaneous reactions (q = 1) of 
two particles this result was shown to be rigorous [15,26]. 

For slow annihilation rates the case k = 2 exhibits a regime of the mean-field like 
behavior (c ~ £ _1 ), followed by a crossover to the fluctuating-dominated asymptotic 
power law (c ~ t -1 / 2 ). This is illustrated by the q = 0.01 data in Figure 2(a). These 
data were averaged over 200 independent Monte Carlo runs for periodic lattices of 
6 x 10 4 sites. Although finite size effects were not severe, rather large lattice sizes were 
necessary to explore large-time behavior while keeping the particle number N(t) large, 
as discussed earlier. The mean-field regime became more pronounced as q decreased. 
For instance, for q = 0.001 (not shown here) it extended over more than three decades 
in the t variable. 

However, the asymptotic behavior of the k = 2 reaction-diffusion system is always 
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fluctuation-dominated, a = 1/2, and the anisotropy introduces no detectable changes in 
the density. Interestingly, the initial density is "forgotten" before the onset of the mean- 
field-regime behavior (and of course the initial short-range correlations are completely 
"washed away" in the fluctuation-dominated asymptotic limit). This property is shared 
by the k = 4 results; see Figure 2(b). 

For k = 4 the asymptotic behavior was found to be mean-field for both isotropic 
and anisotropic hopping. The density variation of the fully anisotropic case (h = 1) was 
in close numerical agreement with the results obtained for isotropic hopping (h = 0.5) 
and therefore it can be described accurately by the mean-field rate-equation calculation 
scheme given in [19]. Since the asymptotic particle concentration decay is now slower 
(c ~ t -1 / 3 ), smaller lattice sizes can be used to investigate the behavior for large times. 
For instance, the q = 0.001 data shown in Figure 2(b) represents average over 200 runs, 
for a periodic lattice of 2 x 10 3 sites. 

4. DISCUSSION 

Generally, our simulations have confirmed the phenomenological considerations re- 
garding the asymptotic particle density, presented in Section 2. However, the observa- 
tion that the data for the isotropic and anisotropic cases are so close numerically for all 
times and g-values has not been explained adequately. 

Our dynamical rules introduce correlations between hopping and reaction, however 
they are more appropriate to describe actual chemical systems. Specifically, the par- 
tial reaction probability rates q < 1 might be interpreted as the result of an effective 
potential that particles must overcome in order to annihilate. Thus, collisions between 
particles (due to diffusion) should promote the particle cluster to go over the reaction 
energy barrier. Although these correlated processes are less well described by mean- 
field calculation schemes, for k > 3 the use of this methodology is still justified within 
the fast-diffusion regime (g<l) where the role of such correlations becomes irrelevant 
both in the isotropic and anisotropic case. It is worth pointing out however, that for 
k = 2 and q = 1 , the decoupling between annihilation and diffusion allows for an exact 
solution of the macroscopic particle concentration which is independent of the hopping 
anisotropy. To the best of our knowledge, no way of solving the case k = 2 is presently 
known either with partial annihilation rates or hopping-reaction correlations. 
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Finally, our numerical results were only for systems with periodic boundary con- 
ditions. It is well established that pair correlations in hard-core (nonreacting) particle 
systems with anisotropic diffusion are extremely sensitive to boundary conditions. It is 
also expected that anisotropic diffusion might introduce significant changes in the form 
of unequal-time correlation functions (not studied in our present simulations). Specifi- 
cally, for k = 2 and q = 1 exact analyses are feasible using fermionic techniques. Thus, 
we hope that the present work will set the stage for further studies of the anisotropic 
multiparticle reactions. 
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FIGURE CAPTIONS 

Figure 1. Snapshot of (a part of) the stochastic evolution for the fully anisotropic 
reaction-diffusion system 4A— nnert. Details of the numerical procedure are de- 
scribed in Section 3; specifically, the case illustrated was for the reaction probability 
q = 0.1. 

Figure 2. Macroscopic density c(t) obtained numerically starting from the random 
particle distribution with the initial concentrations p = 0.25, 0.5 and 0.8. The 
open symbols correspond to the fully anisotropic hopping (h = 1). The small dots 
joined by solid lines denote the isotropic-hopping results (h = 0.5). (a) Data for 
k = 2 and reaction probability q = 0.01, averaged over 200 Monte Carlo runs, for 
a periodic lattice of 6 x 10 4 sites, (b) Results for k = 4, q = 0.001, for 2 x 10 3 
lattice sites, averaged over 200 runs. The dashed lines indicate the asymptotic 
slopes corresponding to the predicted power-law behavior described in the text. 
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